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Abstract 

In the present work, we numerically explore the existence and stability properties 
of different types of configurations of dark-bright solitons, dark-bright soliton pairs 
and pairs of dark-bright and dark solitons in discrete settings, starting from the 
anti-continuum limit. We find that while single discrete dark-bright solitons have 
similar stability properties to discrete dark solitons, their pairs may only be stable 
if the bright components are in phase and arc always unstable if the bright com- 
ponents arc out of phase. Pairs of dark-bright solitons with dark ones have similar 
stability properties as individual dark or dark-bright ones. Lastly, we consider col- 
lisions between dark-bright solitons and between a dark-bright one and a dark one. 
Especially in the latter and in the regime where the underlying lattice structure 
matters, we find a wide range of potential dynamical outcomes depending on the 
initial soliton speed. 

Key words: Solitons, Dark solitons. Discrete Nonlinear Schrodinger equations, 
Manakov model. 



1 Introduction 

The setting of multi-component dispersive systems of nonlinear waves is one 
which has had a significant impact on a variety of areas including the non- 
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linear optics of fibers and crystals [T] and the recently blossiming area of 
Bose-Einstein condensates (BECs) in atomic physics [2|. A fascinating exam- 
ple among the complex nonlinear localized structures that can emerge in such 
settings consists of the so-called symbiotic solitons. These are waveforms that 
would not otherwise exist in unary systems, but are nevertheless supported in 
multi-component ones, precisely because of the inter-species interaction (be 
they atomic species in BEC or different frequencies, polarizations or other op- 
tical "species" in nonlinear optics). One of the principal examples of this type 
consists of the dark-bright (DB) solitons in two-component systems with the 
self-defocusing nonlinearity. The latter class of nonlinearities is well-known to 
support dark solitary wave structures |3^. However, in the presence of non- 
linear interspecies interactions, these dark solitons form a localized potential 
which allows the trapping in the form of bound states of bright-soliton-like, 
density bumps leading to the formation of DB waves. 

Remarkably, such DB structures have been experimentally monitored both in 
nonlinear optics and in atomic physics. They were first experimentally created 
in photorefractive crystals through the pioneering experiments of |4j and subse- 
quently their interactions were partially monitored in [5jp] More recently, the 
realization of multi-component atomic BECs [7f5|^ has prompted a number 
of theoretical investigations of DB waves in the latter setting as well starting 
with the work of [lO], which examined the trapped dynamics of DB soli- 
tons in the presence of the parabohc traps relevant to BECs. Subsequently, 
these structures were also extended to more complex settings such as the 
spinor three-component system, where dark-dark-bright or bright-bright-dark 
solitons could arise [TT]. More recently, the interaction between such dark- 
bright solitons p^llTS] and their potential higher- dimensional generalizations 
into symbiotic, so-called, vortex-bright solitons [14j have been also considered. 
Perhaps more importantly, the experiments of [15] were able to provide the 
first observations of such DB solitons and their interactions e.g. with another 
dark soliton. A realization of multiple dark-bright solitons has also emerged in 
the experiments of [16] and the oscillations of dark-bright solitons in a trap, as 
well as the interaction between two dark-bright solitons were experimentally 
monitored in |17j . 

On the other hand, an area of studies that has also received significant at- 
tention both within the realm of nonlinear optics and within that of atomic 
physics concerns the study of nonlinear dynamical lattices. In these, the "evo- 
lution variable" is continuum, while the "spatial variables" are discrete. One 
of the principal reasons for the growth of this field has been the develop- 
ment of optically induced lattices in photorefractive media, such as Strontium 



^ It should be noted that theoretical proposals of symbiotic DB solitons (although 
in a somewhat different setting than what will be considered below) have existed 
much earlier and at least since the work of [6]. 
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Barium Niobate (SBN) [18], and their experimental realization [T9p0] . As a 
result, a remarkable set of nonlinear waves has been predicted and experi- 
mentally observed; these have now been summarized in a number of reviews 
[2T1I22] . In addition, another direction of nonlinear optics that has grown in 
parallel (and has partially preceded the one above) concerns the intriguing 
interplay of nonlinearity and discrete diffraction emerging in fabricated Al- 
GaAs waveguide arrays |23]. The numerous interesting phenomena observed 
therein including discrete diffraction, Peierls barriers, diffraction management, 
gap solitons, vector solitons, and modulational instabilities among others have 
now been summarized e.g. in [2^21f22j . Lastly, over the last decade, nonlinear 
dynamical lattices have emerged as a theme of interest in BECs since droplets 
of such condensates can be trapped in periodic optical potentials often referred 
to as "optical lattices" (see, e.g., [2S] and references therein). Among the prin- 
cipal byproducts of these studies have been the manifestation of dynamical 
instabilities, Bloch oscillations, Landau-Zener tunneling, and gap solitons as 
has been summarized in the reviews [27|l28] . 

Although, to the best of our knowledge dark-bright solitons have not been 
realized in dynamical lattice settings, all of the relevant ingredients illustrating 
their experimental tractability are available. Namely, recent nonlinear optical 
experiments have considered the realization of dark solitons in the presence 
of self-defocusing nonlinearities of both the Kerr type in AlGaAs waveguide 
arrays [22] and of the saturable type due to the photovoltaic effect in lithium 
niobate arrays [SO]- On the other hand, vector dynamical lattice systems have 
also been recently experimentally monitored; see e.g. [3l] and the discussion 
therein. 

From the above background, it becomes clear that the investigation of DB 
waveforms in discrete two-component self-defocusing systems is a timely and 
relevant theme. Our scope in the present work is to numerically explore this 
topic and the corresponding existence, stability and collision dynamics of DB 
states in nonlinear dynamical lattices starting from a well-established limit, 
namely the anti-continuum one of [32j. The latter has been extremely helpful 
in deriving properties of bright and dark discrete solitons in both focusing and 
defocusing lattices, as can be seen in the detailed discusssion of [33]. In what 
follows, we will see that some of the DB states (such as a single DB soliton or 
pairs of DB solitons with dark ones) draw significant parallels with previously 
studied cases such as single-component discrete dark solitons. On the other 
hand, we will also show that other states such as standing wave pairs of two 
dark-bright solitons have some surprising properties which defy the single- 
component intuition. As will be seen the linear stability of such states will 
be critically affected by the relative phase of the bright components. Lastly, 
we will examine the collisional dynamics of two DB solitons, as well as that 
of a DB state with a dark one. The latter will reveal a significant wealth of 
potential outcomes depending on a single parameter, namely the wave speed. 
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Our presentation will be structured as follows. In section 2, we present the rel- 
evant prototypical model setup in the form of two coupled defocusing discrete 
nonlinear Schrodinger (DNLS) equations. Then, in section 3, we explore the 
existence and stability properties of the DB solitons. In section 4, we consider 
their coUisional dynamics. Finally, in section 5, we summarize our findings and 
present our conclusions. 



2 Model Setup 

We consider the system of two coupled discrete nonlinear Schrodinger (CDNLS) 
equations given by [55] : 

lUn + {9ll\Un? + 9l2\Vn?)Un + CAf/„ = 

iVn + {gi2\Un? + 922\Vn\^)Vn + CAK = 0, (1) 

where A is the discrete Laplacian operator(Af/„ = f/„+i + Un-i — ^Un)- 

Our focus in what follows will be on stationary solutions of the system of 
equations ([T|. To that effect we use the standard ansatz 



Unit) = exp{-iAit)un ; Vn{t) = exp(-zA2t)f n, (2) 

with frequencies (propagation constants in nonlinear optics or chemical po- 
tentials in BECs) Ai and A2, and spatial profiles {un} and Using ^ in 
([T|, we obtain the system of coupled difference equations 

Aiu„ + {gu\un\^ + gi2\Vn\'^)Un + CAu^ = 

MVn + {gi2\Vn\^ + g22\Vn\'^)Vn + CAVn = 0. (3) 

The anti-continuum (AC) limit of [32] (see also [33]) concerns the case with 
C = 0. In the absence of the potential, the system of Eq. ^ results into a 
pair of analytically solvable algebraic equations for each site. In what follows, 
and given the nature of the DB solitons, we will restrict our attention to cases 
where the first (dark) component is excited and the second (bright) one is not 
i.e., Un = ±^J—Ai/gu, Vn = (at the tails of the DB soliton) or to ones where 
the first component is inert and the second one is excited i.e., Un = and 
Vn = -^\J —A2/ g22- (at the center of the DB sohton). Once an exact solution 
of the DB type has been constructed at the AC limit, it can be continued 
(numerically, as well as analytically) for finite values of the coupling C. It 
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is worth remarking that throughout the text, unless otherwise specified, free 
ends boundary conditions have been used. 

A relevant question concerning the solutions at finite values of C is that of 
dynamical stability. To that effect, we use linear stability analysis to provide 
a spectral response to this question. The stability is determined in a frame 
rotating with frequency Ai for Unit) and A2 for i.e., we suppose that 

m 

Unit) = exp(-zAit)K + Vnit) = exp(-zA2t)K + ^l^\t)]. (4) 

The small perturbations ^nK'^)^ with k = 1,2, can be expressed as 

Cn\t) = a„exp(ia;t)+6„exp(-iw*t), Cn\t) = c„ exp(2wt)+(i„ exp(-iw*t), 

(5) 

leading to the linear stability equations 
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with 

i^l,n = Ai + 2gii\Un\^ + gi2\vn\^ + A, 
^2,n = A2 + 2g22\Vn\^ + gi2\Un? + A. 

More specifically, it can be proved [33] that a solution is spectrally stable if 
all the u eigenvalues (the so-called eigenfrequencies of the system) are real 
numbers. Two type of instabilities are observed in this system: (i) exponential 
instabilities, represented by a eigenfrequency (w) pair with zero real part and 
non-vanishing imaginary part and (ii) oscillatory instabilities, accounted for 
by a Hamiltonian Hopf bifurcation, and represented by quartets of eigenfre- 
quencies with non-zero real and imaginary parts. 

Throughout this paper we have assumed gu = —0.97, (722 = —1-03 and gi2 = 
— 1, chiefiy motivated by the experimental discussed interspecies interaction 
ratios of fTl'S]- AH the calculations have been performed assuming (without 
loss of generality) Ai = 1, and for the frequency A2 > we have explored the 
range A2 G [0, 1.5], as representative of the different parametric regimes. 
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3 Stability properties of dark-bright multi-soliton solutions 

3.1 Solutions at the anti- continuous limit 

In this subsection, we examine a representative set of structures which are 
initialized at C = 0. In the AC limit, we use the following convenient notation 
to represent our solutions: by indicating = ±1, we denote an excited site 
of Un = i^— Ai/^fii, while similarly for the second component {;„ = ±1 is 

used to denote f„ = iy^—A^T^- In this notation, we have considered the 
following (AC-limit) solutions as seed points for our numerical continuations: 

• Dark-bright one-soliton solutions (DB 1-S). 

Un = sgn(n), Vn = Sn,o, (7) 

where sgn(n) denotes the sign of n and it is implied that the solution 
assumes the value Uq = 0. 

• Dark-bright two-soliton solutions (DB 2-S). We have considered, as a rep- 
resentative case, that where the distance between the excited bright spots 
is of 2 sites: 

■Un = 1 - S\n\,l - 2(5n,0- (8) 

In this setting, we distinguish two configurations for the bright part: 
the in-phase bright case with 

Vn = Sn-l + Sn,l, (9) 

and the out-of-phase (anti-phase) case with 

Vn = Sn-1 - Sn,l- (10) 

• Dark-bright and dark two-soliton solutions (DB+D). Again, we have con- 
sidered that there are two sites separating the vanishing sites of the m„ field, 
one of which is associated with the dark-bright soliton, while the other is 
associated with the dark one. 

Un = 1- 5\n\,l - 25nfi, Vn = Sn-1 (H) 

Figure [T] shows some examples of the profile of the above mentioned configu- 
rations. 

It is interesting to note that at the AC limit, not only the existence theory 
is transparent but also the nature of the linearization spectrum. More specifi- 
cally, notice that all the DB states involve configurations for which UnVn = 0, 
hence the above 4N x 4N (where N is the number of lattice sites considered) 
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Fig. 1. Profile of the configurations analyzed in the paper. In all cases, C = 0.02, 
Ai = 1 and A2 = 0.15. 

matrix splits into two diagonal 2N x 2N blocks for each of the components. 
The background states (away from the DB center) provide each of these com- 
ponents with a spectral dispersion relation of the form: 



u = ±^4Csin2(|) [4Csm\^) + 2Aij (12) 

between the linearization frequencies u and wavenumbers q; this relation stems 
from the dark (first) component. On the other hand, the bright (second) com- 
ponent similarly yields a dispersion relation for the background state (away 
from the center of the DB) 



u; = ± (^A, _ A2 + 4Csin2(^)^ (13) 

For the excited sites of the dark-bright soliton (where the first component is 
vanishing i.e., dark and the second component is non- vanishing i.e., bright), 
we obtain a pair of eigenfrequencies satisfying: 



c^ = ±f-— A2 + A1) (14) 

V 922 J 

from the dark component and a pair of zeros from the bright one. The funda- 
mental stability question concerns, of course, the fate of these eigenfrequencies 
when C 7^ and is examined systematically in our numerical computations 
within the following subsection. 
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Fig. 2. Existence and stability properties of the dark-bright one-sohton solutions as 
a function of coupling C and frequency A2, for the case of Ai = 1. 

3.2, Numerical continuations for arbitrary coupling 

Having considered the existence and stability properties of the DB waveforms 
at the AC hmit of C = 0, we now turn to the case of finite couphng C. In what 
follows, we examine for the configurations discussed in the previous subsection, 
the findings pertaining to their existence when a continuation from C = is 
performed, calculating for each solution the eigenfrequencies stemming from 
the relevant stability problem ([6]). In every case, Ai = 1 is fixed whereas A2 
is chosen in the interval [0, 1.5]. 

3.2.1 Dark-bright one-soliton solutions 

Three main behaviours are observed separated by the critical values A2 = 0.51 
and A2 = 1. An overall view for the existence and stability properties of these 
solutions, as a function of coupling C and frequency A2, is illustrated in Fig. |2] 

For A2 < 0.51, stable solutions appear only within a small parametric interval 
near the AC limit. In this case, past a (A2 dependent) critical point an insta- 
bility arises due to the collision of a localized eigenmode with the top of the 
dark component linear mode band, which leads to a Hopf bifurcation. Addi- 
tionally, the amplitude of the bright spot decreases as C is increased (for fixed 
A2) vanishing above a critical value of C. For small values of A2, the solution 
vanishing takes place before the destabilization, implying that the DB soliton 
is stable in the whole range of C in which it exists with a nontrivial bright 
component. Fig. |3] (top) shows for a prototypical value of A2 within this pa- 
rameter range, the dependence of the eigenmode spectrum with respect to C. 
The figure also features the comparison of the relevant eigenfrequencies with 
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Fig. 3. Real part of the relevant spectral eigenfrequencies (left panels) and imaginary 
part of the eigenfrequency spectrum (right panels) versus C of the DB one-soliton 
solutions. In both cases, Ai = 1; the top panels are for A2 = 0.15, while the bot- 
tom ones for A2 = 1. The dashed lines in the left panels correspond to analytical 



predictions (12)-(14). Recall that instabilities emerge when Im(a;) 7^ 0. 



the analytical predictions 
C. Recall that per the ab 
(pertaining to the "dark 
with the expanding band 
fact leads to the observed 
critical point in C. When 
the bottom panels of Fig. 
to 0. 



(12)- (14) showing an excellent agreement for small 



ove analytical predictions there is an internal mode 
part" of the DB object) which collides at finite C 
of phonon spectrum growing away from zero. This 
instability of the DB solitary wave past the relevant 
A2 is increased, as it is evident from both Fig. |2] and 
[3} the critical C for the instability emergence tends 



For A2 G (0.51, 1] solutions exist for every value of C. Now, once the lo- 
calized mode collides with the linear mode band, new parametric regions of 
stabilization appears as "bubbles" caused by a Hopf bifurcation cascade, as 
shown in Fig. |4j The size of the instability bubbles decreases and the num- 
ber of them increases as the number of particles increases. This is the so 
called boundary-induced stabilization and is a finite size effect discussed e.g. 
in [35] . Figure |5] illustrates, as an example, two stability diagrams obtained 
with N = 101 and N = 10001. The instability bubbles' size is almost inde- 
pendent of the boundary conditions used, whereas the number of bubbles is 
independent of them, as shown in Fig. |6} Fig. [3] (bottom) shows the eigen- 
modes spectrum for an example within this range. This behaviour resembles 
to that of a dark soliton |35j|. Contrary to the above mentioned case, in this 
parametric regime, there is no vanishing of the bright spot and genuine DB 
solutions exist for every value of C. Finally, for A2 > 1, the DB soliton is 
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Fig. 4. Real part of eigenfrequencies versus coupling C for dark-bright one-soliton 
solutions with A2 = Ai = 1. The right panel is a zoom of the left panel where the 
the existence of a Hopf bifurcation cascade is highlighted. 
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Fig. 5. The imaginary part of the eigenfrequencies is plotted versus the coupling C 
for A2 = 0.9. The dependence on the system size is illustrated: N = 101 is chosen 
in the top panel whereas = 10001 is taken in the other ones. The bottom panel 
highlights a zoom-in of the middle panel into a parametric range where bubbles 
appear. 

unstable for every value of C. 



3.2.2 Dark-bright two-soliton solutions 

Here we consider both the in-phase and anti-phase cases. 

In the in-phase case, we observe that the DB solitons behave structurally 
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Fig. 6. The imaginary part of the eigenfrequencies is plotted versus the couphng 
C for A2 = 0.9 and N = 101. The (in) dependence of the stabihty bubbles on the 
boundary conditions is highlighted. The top panel corresponds to free end boundary 
conditions, the middle panel to anti-periodic boundary conditions and fixed ends are 
shown in the bottom panel. Notice that the imaginary parts of the eigenvalues are 
almost identical for fixed ends and anti-periodic boundary conditions (the differences 
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Fig. 7. Real part of the relevant eigenfrequencies (left panels) and imaginary part of 
the full spectrum (right panels) versus C of the DB two-soliton solutions. Dashed 



lines in left panels correspond to analytical predictions (12)-(14) 



in a similar way to the case of the single dark-bright soliton solutions with 
A2 < 0.51 (left region of Fig. [2]), that is, for small C the solutions are sta- 
ble up to a critical value where they destabilize, with their bright component 
disappearing above a second critical value of C; additionally, stability bub- 
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Fig. 8. Existence and stability properties of the dark-bright two-sohton solutions as 
a function of coupling C and frequency A2, with the fixed value Ai = 1. The top 
panel shows the in-phase case, while the bottom one illustrates the anti-phase case. 

bles do not exist. Similarly, there do not exist regions where the DB pair 
solution persists for arbitrary C . The latter is natural to expect since the in- 
phase structure induces repulsion among the bright component^^l this fact, 
combined with the repulsion among the dark components, does not allow the 
possibility of sustaining stationary such solutions near the continuum limit. 
Moreover, we observe that the maximum value of C does not depends mono- 
tonically with A2, having a maximum around A2 = 0.42. For A2 below this 
critical value, and upon increase of C, the bright spot amplitude vanishes 
similar to the DB 1-S case but now the solution that exists above the critical 
value of C consists of 2 dark soliton waveforms; however, when A2 is above 
the critical value, the continuation fails past a certain critical C. This is due 
to the existence of an inverse Hopf bifurcation, which may be caused by the 
interaction between the two peaks of the bright component, similarly to what 
is observed in 1361. 



On the contrary, the out-of-phase case solutions are always exponentially un- 
stable. Additionally, we observe a monotonically increasing behaviour of the 
maximum value of C with A2 up to A2 ~ 0.65. In this parametric range, 
the bright spot vanishes above the critical coupling. On the other hand, for 
0.65 ^ A2 < 0.75, the curve is almost horizontal and the solutions experience 
an inverse Hopf bifurcation. Finally, DB 2-S exist for any value of C whenever 



^ This is an important structural characteristic of the self-defocusing nature of 
our nonlinearity. Namely, contrary to what is standardly known for the case of 
self- focusing/attractive interactions [1], bright solitons here repel when they are in 
phase, but attract when they are out of phase; cf. also with the continuum case 
discussion of Ref . [13] . 
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Fig. 9. Existence and stability properties of the dark-bright and dark two-sohton 
solutions as a function of coupling C and frequency A2, with the fixed value Ai = 1. 

A2 > 0.75. These findings are summarized in Figs. [7|and|8j Let us mention in 
passing that somewhat surprisingly, the stability for small C is the opposite 
than for the case of bright solitons which are separated by an even number 
of sites (i.e. in the latter case, the bright 2-solitons in-phase are unstable for 
any coupling whereas the bright 2-solitons in anti-phase are stable for small 
coupling) [Sn] ; for a corresponding discussion in the multi-soliton interactions 
in the 2d case, see e.g. [57] . 

3.2.3 Dark-bright and dark two-soliton solutions 

The dark-bright and dark two-soliton solutions behave structurally in a similar 
way to the DB 1-S or DB 2-S in-phase solutions. There are however some 
important differences. The maximum value of C (as a function of A2) for 
which such solutions exist has a maximum for A2 = 0.70. Below this value, the 
amplitude of f 1 < (i.e. the amplitude of the absent bright spot), increases its 
absolute value, and, above the maximum C the DB+D soliton transform into 
a DB 2-soliton solution in anti-phase. For A2 > 0.70, the soliton experiences 
an inverse Hopf bifurcation at the critical C. Figure |9] illustrates the existence 
and stability properties of the dark-bright and dark two-soliton solutions as a 
function of the coupling C and the frequency A2. 

3.2.4 Summary of existence and stability results 

The previous findings for dark-bright two-soliton solutions and for the dark- 
bright and dark two-soliton pair solutions were obtained for a separation of 
two sites between holes or peaks. However, the results are similar when the 
separation increases up to 3 sites. The relevant effects on the lattice decay 
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exponentially with the separation, so that for large values of this separation, 
the DB structures behave as isolated ones and the stability properties are 
similar to the dark-bright one-soliton solutions. 

Arguably, the most interesting of the above results concerns the stability of 
the DB two-soliton state and especially the independence of this result on the 
distance between the DB waves; instead, the result seems to solely (and criti- 
cally) depend on the relative phase of the bright components. It is relevant to 
contra-distinct this to the case of simply two-dark solitons of the form defined 
above which would be stable (at least close to the AC limit) independently 
on the distance d, and also to the case of two bright gap solitons, considered 
before in |38] (see also [33]) whose stability would depend on both the distance 
d and the relative phase. In fact, in the latter case of two bright gap solitons 
for defocusing nonlinearities, using the so-called staggering transformation to 
convert the defocusing dynamics into focusing one, it is straightforward to 
establish [3H] that in-phase excited sites, separated by an even number of sites 
are unstable, while by an odd number of sites are stable (for small C), while the 
situation is reversed for anti-phase excitations (stable for an even site distance 
and unstable for an odd one). On the other hand, here we observe that the 
"dressing" of the interaction of the bright components via their dark partners 
leads to stability for in-phase excitations and instability for out-of-phase ones. 
At a qualitative level, one can present the following straightforward argument 
towards understanding this effect: it is well-known (see also the discussion 
below regarding DB soliton collisions) that DB solitons with in-phase bright 
components repel each other while DB states with out-of-phase bright exci- 
tations attract each other; see e.g. the fundamental experimental work of [5], 
and the very recent theoretical discussion of [39]. As a result the state with 
out-of-phase bright excitation will be structurally less robust e.g. if the DB 
pair is moved closer, the attraction becomes stronger and the tendency is to 
deviate further from the equilibrium. If the pair is moved further away, the at- 
traction becomes weaker, and again the further deviation from equilibrium is 
favored. The contrary is true for the in-phase excitation. The repulsion of the 
bright components, if they are brought closer or drift away from the equilib- 
rium, favors in both cases their return to equilibrium indicating its structural 
robustness. 



4 Collisions 

Having examined the stability of the DB soliton configurations, we now turn 
to a more systematic examination of collision events involving two DB waves 
or a dark-bright and a dark one. All the numerical computations below have 
been performed perturbing only the bright parts of the corresponding solu- 
tions. In the case of dark-bright and dark two-soliton pairs, the perturbation 
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only affects the single bright spot, and for the other cases the perturbations 
affect symmetrically each bright wavepacket so that they can move in oppo- 
site directions. For most of the collisions we have taken N = 402, Ai = 1, 
A2 = 0.7 and for each value of C, the values of a (the parameter controlling 
the "kick" oc e*"("-"«), where no is the center of the soliton), have been taken 
in the interval [0.01,0.30]. The velocity of the dark-bright soliton increases 
when a increases. The solitons are initially separated by a distance of 200 
lattice sites. Notice that only collisions at intermediate and high values of C 
have been considered, as for low values of C the solitons are not mobile due 
to discreteness effects which impede their traveling. We have observed that, 
for A2 = 0.7, which is the case analyzed in most of the cases below, C must 
be higher than 1.2 if a = 0.01 and 0.9 if a = 0.3. The critical coupling needed 
for mobility decreases with A2 and a. 



4-1 Collisions between a dark-bright and a dark soliton 

(1) Collisions at intermediate coupling 

First of all, we consider collisions for C = 1.23. The possible outcomes 



of such events are illustrated in Figs. 10 11, obtained with increasing 
values of a. This value of C has been chosen so that it is "intermedi- 
ate" between the small coupling regime (of low soliton mobility) and the 
high coupling regime where the collisions are essentially tantamount to 
the ones in the continuum model (and discreteness is irrelevant). In this 
intermediate range, the discreteness plays a crucial role in determining 
the coUisional outcome and provides a rich variety of potential scenarios. 
The collision events are illustrated by means of space-time diagrams rep- 
resenting the contour plots of the densities of the two fields, {|ttnP} and 
{|w„p}, respectively. 

For a small enough the static dark remains at rest after the collision. 



and the dark-bright is reflected from it, as shown in Fig. 10 for a = 0.08. 
For larger values of a such as 0.12 of the middle panel of Fig. 10 , the dark- 
bright may rebound but also concurrently leads to the motion alongside it 
of the stationary dark soliton, formulating an interesting co-propagating 
dark-bright and dark soliton pair. This apparent effect of discreteness 
is not only absent for higher couplings C (as we will see below), but 
also is especially unexpected given the apparently repulsive nature of 
the interaction between the dark-bright and the dark soliton for large 
C (i.e., in the continuum limit), discussed below. This scenario appears 
for a small interval of a values. The explanation that we offer for this 
effect is the following: for the present value of C and the relevant lattice 
size = 402, the dark soliton is, in fact, dynamically unstable. While 
this instability does not manifest itself until the collision event, the latter 
appears to trigger the instability by providing a significant perturbation 
to the stationary dark soliton. As a result, after the collision the dark 
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soliton becomes subject to the oscillatory instability, in turn giving rise 
(as observed previously) to its motion; see also pU] . 

For increasing values of a the static dark always gets into motion but 
some different scenarios appear: 

(a) The dark-bright solitary wave rebounds launching out away the static 
dark, as illustrated for a = 0.13 in Fig. [lOj 

(b) The dark-bright wave may remain nearly stopped, as for a = 0.24 in 



the left panel of Fig. 11 



The dark-bright may continue traveling along the incident direction, 
as the dark (initially stationary) soliton does; an example of this for 
a = 0.245 can be found in the middle panel of Fig. 11 



(d) Lastly, we observe another unusual dynamical scenario unfold in the 



right panel of Fig. 11 here, apparently the dynamical instability of 
the stationary dark soliton manifests itself prior to the arrival of 
the DB soliton (and the ensuing collision event). As a result, after 
the collision the dark-bright continues traveling without essentially 
changing its velocity while the dynamically unstable dark soliton 
continues traveling in the opposite direction. 



Dark fa=0.06) 



□ark[a=0.135 




„ 5200 



Bright [n=0.12J 



Bright [n=0.12) 



Bright [a=0.13) 




Fig. 10. Three dark-bright+dark 2-S collisions corresponding to increasing values of 
a; Q = 0.08 on the left, a = 0.12 in the middle and a = 0.13 on the right. Relevant 
parameters are chosen as: C = 1.23; Ai = 1; A2 = 0.7. 

(2) Collisions at higher coupling 

Some prototypical examples of collisions at higher coupling are also 



given for C = 3 in Fig. 12 These are essentially for a regime where dis- 
creteness is irrelevant and the dynamics chiefly represents the near-elastic 
collisions of the continuum limit (which, however, can still be fairly com- 
plex, exhibiting beating and breathing phenomena as illustrated in j41]). 
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Fig. 11. Same as Fig 10, but for larger values of the initial kick a (0.24, 0.245 and 
0.25, respectively on the left, middle and right), corresponding to larger speeds. 

Nevertheless, we show them here to showcase the repulsive nature of the 
interaction between the dark-bright and the dark soliton. It is clear that 
the low velocity collision of the left panel of the figure signals a repulsive 
occurrence whereby the centers of the two solitons never "touch" . There 
is a critical kinetic energy which enables the DB soliton to balance the 
repulsion induced from the stationary dark soliton. Close to this critical 
point, there is a prolonged merger of the two waves, prior to their separa- 
tion as evident in the middle panel of the figure. Finally, the right panel 
shows a supercritical case, where the DB soliton possesses enough kinetic 
energy to overcome the repulsive barrier posed by the dark stationary 
state. In that case, the dark remains stationary, while the dark-bright 
passes through it. It should be noted that in these events, the large value 
of C mitigates the dynamical instability (which now has a very small 
growth rate as the continuum limit is approached), hence some of the 



more elaborate occurrences of Figs. 10 11 are absent here. 



4-2 Two dark bright solitons collisions 



In order to observe the collision scenario between two DB solitons, we fix 
C = 1.55. This value of C has been chosen so that it is "intermediate", as in 
the above examined collisions between a dark-bright and a dark soliton. 

First of all, we consider slow soliton collisions. The parameter a should be 
small enough so that the kinetic energy does not overwhelm the details of the 
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Carh (tt=0.03) 



□ark (0=0.1 4.) 



□ark(a=0.141) 




collisions. For these low velocities dark-bright solitons with in-phase bright 
components repel each other while dark-bright solitons with out-of-phase bright 
excitations attract each other. These results agree with the experimental work 
of [5] and the theoretical work of [13] which discusses the nature of the inter- 
action. 



C=1.55 



o" Dark (a=O.01) in-phasB ^^u' Dark ([i=0.[)1) anti-phase 




-40 -33 20 40 -40 -20 20 40 



Q^i Bright (a=0. 01) in-phasg ^^^j" Brigtnt (a=0. 01 ) anti-phase 




Fig. 13. An example of dark-bright two-soliton collisions. (Left): in-phase case. 
(Right): anti-phase case. Parameters: C = 1.55; Ai = 1; A2 = 0.7 and a = 0.01. 
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As an example, taking a = 0.01, Fig. [13] illustrates a dark-bright two-soliton 
collision event for the in-phase case (left column) and for the anti-phase case 
(right column), with A2 = 0.7. We observe repulsion of the dark components 
when the bright components are in phase and attraction when they are out of 
phase. 



If the velocity is increased, the following modifications are observed. On the 
one hand, in the in-phase bright component case, the repulsion of both the 
dark and bright components is mitigated by the increase of kinetic energy 
(i.e., the turning points of the relative motion occur closer to n = 0); on 
the other hand, for the anti-phase bright component, the kinetic energy and 



attractive interaction cooperate in bringing the solitons together. Fig. [14] shows 
an example of the outcome of both cases for a = 0.3. 



□ark 1(1=0.3] in-phasfl 



□ark (a=0.3) anti-phi 





Bright (ft=0.3) in-phas9 



BriginT (a=0.3) anti-phase 





Fig. 14. Same as Fig. 13, but for a = 0.3 that corresponds to a larger speed. 



In order to discern if the dark components of the two DBs with vr out of phase 



bright components cross each other or not, we have traced in Fig. 15 the dark 
component density profile |MnP at the time where both colliding solitons are at 



the minimum distance for the right panels of Figs. 13 and 14 We can observe 
the existence of two dips instead of a single one, which is a clear indication 
that the dark components of the solitons do not cross each other. In fact, as 
discussed also in |12], it is only when the two dark soliton centers coincide 
that they pass through each other (i.e., that they overcome the finite barrier 
of repulsion between them). 

Ref. [13j considers the interaction between two DBs with vr out of phase bright 
components as the norm of the bright component is modified. There, the fading 
attraction effect between the bright components (as their norm decreases) 
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Fig. 15. Density profile of the dark components of the colhding sohtons of right 



panels of Fig. 13 (left) and Fig. |14| when the distance between them is the shortest. 



enables the emergence of a repulsive character of the interaction due to the 
dark components. In order to check if this phenomenon also occurs in our 
system, we have performed some simulations varying A2 (increasing A2 is 
equivalent to increase the norm of the bright component) showing that there is 
a critical value of A2 below (above) which the dark components repel (attract) 



each other. This phenomenon is illustrated in Fig. 16 



Finally, we have checked that, contrary to the dark-bright and dark solitary 
wave collisions, no fundamental qualitative differences are observed when a 
high coupling in implemented. That is, collisions at intermediate and strong 
couplings C roughly exhibit the same principal phenomenology (hence, we do 
not show any of the latter here). 



5 Conclusions and Future Challenges 



In the present work, we offered a numerical perspective on the existence, stabil- 
ity, dynamics and interactions of single and multiple dark-bright solitons not 
only between them, but also with dark solitons in the setting of nonlinear dy- 
namical lattices. The presence of discreteness enabled a number of interesting 
variants to the previously studied continuous setting. Some of these variants 
could be anticipated on the basis of analogous studies of dark solitons [33]. E.g., 
the dark-bright single solitons could be stable only up to a critical coupling 
(although finite-size effects could provide additional stabilization). The same 
thing would be true for pairs of dark-bright and dark solitons. The collisions 
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Fig. 16. Comparison between 2 DB soliton collisions when A2 is varied (A2 = 0.7 
on the left, while A2 = 0.52 on the right; C = 5 and a = 0.1 in both cases). It is 
observed that below a critical value, the dark component repulsion can be deemed 
to be stronger than the anti-phase bright component attraction and hence creates 
a potential barrier that, in turn, induces a finite separation "turning point" in the 
interaction between the solitary waves. 

of the latter pairs indicate the repulsive interaction between the dark-bright 
and the dark states. On the other hand, a number of interesting novel fea- 
tures were revealed. Relevant examples include the generic instability of the 
dark-bright for a sufficiently strong bright component, the generic instability 
of out-of-phase two-DB soliton pairs, and the generic stability (again for large 
enough bright component) of in-phase two-DB soliton pairs. Discreteness was 
also found to play a significant role in affecting collisions for "intermediate" 
coupUng strengths (not too small, so that the excitations would get pinned, 
nor too large so that they would interact in a quasi-continuum way). There, its 
potential combination with the dynamical instability e.g. of stationary dark 
solitons could yield a variety of interesting coUisional outcomes, especially in 
the collisions of dark-bright solitary waves with dark ones. 



There are numerous directions in which one can consider expanding the con- 
siderations presented herein. One possible example is an analytical attempt 
to calculate the linearization eigenvalues of two-DB soliton pairs, to explicitly 
prove the stability conclusions numerically inferred herein (and to obtain some 
mathematical intuition on these findings). Another possibility from a numeri- 
cal viewpoint is to extend the considerations presented herein to vortex-bright 
solitons (or to dark-bright rings) in higher-dimensional settings, by analogy 
to the continuum form of such symbiotic entities recently considered in [13]. 
There one can potentially envision situations where the symbiotic structure 



21 



may be stable although neither of its constituents is robust by itself. Potential 
experimental realizations of dark-bright solitons and soliton pairs in waveguide 
arrays (such as lithium niobate ones) could be an excellent source of testing 
of the relevant ideas in realistic settings. 
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